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Abstract. The external electric field deforms flaccid phospholipid vesicles into 
spheroidal bodies, with the rotational axis aligned with its direction. Deformation 
is frequency dependent: in the low frequency range (~ 1 kHz), the deformation is 
typically prolate, while increasing the frequency to the 10 kHz range changes the 
deformation to oblate. We attempt to explain this behaviour with a theoretical model, 
based on the minimization of the total free energy of the vesicle. The energy terms 
taken into account include the membrane bending energy and the energy of the electric 
field. The latter is calculated from the electric field via the Maxwell stress tensor, 
where the membrane is modelled as anisotropic lossy dielectric. Vesicle deformation 
in response to varying frequency is calculated numerically. Using a series expansion, 
we also derive a simplified expression for the deformation, which retains the frequency 
dependence of the exact expression and may provide a better substitute for the series 
expansion used by Winterhalter and Hclfrich, which was found to be valid only in 
the limit of low frequencies. The model with the anisotropic membrane permittivity 
imposes two constraints on the values of material constants: tangential component 
of dielectric permittivity tensor of the phospholipid membrane must exceed its radial 
component by approximately a factor of 3; and the membrane conductivity has to be 
relatively high, approximately one tenth of the conductivity of the external aqueous 
medium. 
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1. Introduction 

Rapid development in biotechnology induced a growing interest in the influence of the 
AC electric field on biological cells, phospholipid vesicles and colloidal particles pQ , which 
instigated theoretical investigations in this area [2j [3]. Among the systems studied, 
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phospholipid vesicles are often chosen as the model system of choice due to their well- 
defined structure. 

In the late 1950s, Schwan [I] began a series of studies of the influence of the electric 
field on biological cells, modelling them as simple geometric shells with given electric 
properties. In the early 1970s, Helfrich included the then newly developed elastic theory 
for lipid bilayers j5] into the treatment of the effect of the external electric field on 
phospholipid vesicle shape [6J. The first experiments with phospholipid vesicles in an 
electric field which followed several years later [7j seemed to confirm the claim from 
the previous theoretical analysis |j: a 2 kHz external AC electric field deforms flaccid 
phospholipid vesicles into prolate spheroids. Winterhalter and Helfrich [8] extended the 
earlier theoretical treatment [Bj; their model allows for a finite electrical resistance of 
the bilayer, an AC electric field, and includes Maxwell stresses inside the membrane. 
The authors treated the spherical vesicle as a lossy dielectric immersed into a medium 
which was also treated as a lossy dielectric. They examined the case a m /a w <C e m /e w , 
where a m , a w , e m and e w are the conductivities and permittivities of the membrane 
and the aqueous medium, respectively. Employing series expansions, they obtained a 
simple expression for vesicle deformation and predicted a conducting regime of vesicle 
behaviour at low frequencies (u <C o- m /e m ), where the field does not penetrate the 
vesicle interior, a dielectric regime at higher frequencies (u ^> cr w /e w ), where the field 
penetrates the vesicle interior, and an intermediate regime in between the two. In all 
three regimes, they obtained a prolate vesicle shape. 

In contrast to that, experiments with varying frequency P, [10] have demonstrated 
that a vesicle undergoes a prolate-to- oblate shape change when the frequency of the 
applied field is increased from the 1 kHz range into the 10-kHz range (figure Hk). 
Figure [lb shows the dependence of the ratio of vesicle semiaxes c/a for the shown 
vesicle on the frequency of the applied electric field, with c being the semiaxis parallel 
to the field, and a the semiaxis perpendicular to it. 

In this paper, we take the work of Winterhalter and Helfrich as our starting point 
and proceed to extend it by allowing for an anisotropy of the permittivity of the 
phospholipid bilayer. Since the model of Winterhalter and Helfrich departs from the 
experimental data in the high-frequency region by failing to predict oblate shapes, we 
sought a modification that would affect its behaviour in the high-frequency range, where 
the behaviour of the system is governed solely by the permittivities. 

Recently, several papers have been published which take into account the anisotropy 
of the membrane dielectric permittivity and the membrane electric conductivity, which 
arise due to adsorbed hydrophobic ions on the membrane [HJ H2] . Ambjornsson and 
Mukhopadhyay [13] have provided a solution for the electric potential for a general 
ellipsoid coated with a dielectrically anisotropic coating; the same group also examined 
a possible mechanism for the dielectric anisotropy in the membrane at frequencies where 
molecular resonances are important [14]. Simeonova and Gimsa [15] have extended 
the description of the phospholipid vesicle membrane to a three-layer shell to take into 
account the fact that the dielectric anisotropy occurs only in the phospholipid headgroup 
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Figure 1. (a) Giant phospholipid vesicle is deformed into a prolate shape (left) at an 
applied AC field of 1 kHz, and into an oblate shape at 10 kHz (right). The conductivity 
of the aqueous medium was 1.3 U-S/cm. An electric field of ~ 2 ■ 10 4 V/m was applied 
in the horizontal direction. The bar represents 10 |a.m. (b) Deformation of the same 
vesicle, expressed as the value of the semiaxis c (aligned with the direction of the electric 
field) as a function of the circular frequency u> of the external AC electric field. The 
dashed line corresponds to the value at which the vesicle is spherical; greater values of 
c correspond to prolate shapes, smaller to oblate. The dotted line is discussed further 
in the text. The three consecutive runs yield 2394±75, 2427±75 and 2371±87 urn 3 , as 
estimates for the vesicle volume, respectively, which indicates that the vesicle volume 
did not change significantly during the experiment. Details on experimental setup are 
given in [20] , 

layer, while the middle layer of hydrocarbon chains is mainly isotropic. It also needs to 
be noted that transitions between prolate and oblate shapes in response to a frequency 
variation of the applied electric field have been examined theoretically for cases where 
the electrical properties of the aqueous solution inside the vesicle differ from those in 
the vesicle exterior pT6], ITT], and recently, an experimental morphological phase diagram 
has also been obtained |18j . A review of the studies of the influence of the electric field 
on phospholipid vesicles was recently summarized in a paper by Dimova and coworkers 

While the above papers provide an analytical expression for the electric potential 
[13] or a thorough numerical analysis of the Clausius-Mossotti factor [TH [T2J, [15], they 
do not focus on an analysis of parameters at which either prolate or oblate deformation 
is possible in the high- and low- frequency limit, which is the aim of this paper. 

The paper is organized as follows. In section 2 we first develop an extension of the 
model of Winterhalter and Helfrich by taking into account the anisotropy of membrane 
permittivity. In section 3 we show that the model presented here exhibits frequency- 
dependent prolate-to-oblate shape transition of a phospholipid vesicle similar to the one 
observed in the experiment. In the expansion valid for large vesicles, we derive the 
neccessary conditions for a prolate or oblate deformation in both the high- and low- 
frequency limit, and analyse the dependence of the prolate-to-oblate shape transition 
frequency on material constants. Finally, in section 4, we discuss some of the issues 
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which arise from the results. 
2. Theoretical framework 

The equilibrium vesicle shape in an external electric field is calculated as the shape with 
the minimal total free energy, consisting of the vesicle bending energy and the energy 
due to the electric field. The calculation was conducted for a vesicle with constant 
volume, while the necessary area increase for the deformation stems from the flattening 
of thermal fluctuations of the vesicle by the electric field [211 [22] . The analytical 
expression for the membrane bending energy is well known [5], and the change of the 
energy due to the electric field is calculated as the work done by the force of the electric 
field while deforming the vesicle [8]. In order to evaluate it, we first have to calculate 
the electric field. The present treatment is limited to small deviations of vesicle shape 
from the sphere, thus the electric field in the presence of a spherical shell is computed. 
Both the aqueous solution inside and outside the vesicle and the vesicle membrane are 
treated as lossy dielectrics. The aqueous solutions inside and outside the vesicle have 
identical electrical properties. Dielectric permittivity of the membrane is treated as 
anisotropic, with the component aligned with the normal to the membrane differing 
from the components perpendicular to this direction. Both the dielectric permittivity 
and the electric conductivity of the aqueous solution as well as the electric conductivity 
of the membrane are treated as scalars. 

The applied electric field E introduces a single distinct axis into the system, so the 
treatment can be limited to axially symmetric shapes. This eliminates the dependence 
on the longitude angle <ft, if the polar (z) axis is chosen parallel to the applied field. 
Thus, the vesicle surface can be parametrized as 

r(9) = s + s(9) , (1) 

where \s(9)\ <C s , s denoting the deformation independent of the polar angle 9. In an 
absence of deformation (|s(0)| =0), so equals the radius of the undeformed sphere r . 
The deformation is independent of the sign of the electric field, thus it is proportional 
to E 2 in the lowest order. As the field itself is proportional to cos 9, one can expect a 
deformation coupled with a field to be proportional to cos 2 9. In terms of expansion into 
spherical harmonics, this limits us to a sum of even terms. Retaining only the terms 
proportional to E 2 or lower, a quadrupolar term remains, where s(9) equals the second 
Legendre polynomial: s(9) = |s 2 (3cos 2 6' — 1), with s 2 being a measure for the extent 
of deformation. Positive values of s 2 indicate prolate deformation, negative oblate. 

If the membrane area is to be locally conserved, a quadrupolar displacement 5r r in 
a radial direction must be accompanied by a tangential displacement Sr g [S] : 

Sr r = ^(3cos 2 fl- l)s 2 (2a) 
5r e = — cos 9 sin 9 s 2 ■ (26) 

A requirement for a local area conservation assures that the membrane stretching is 
independent of polar angle 9. The total membrane area expansion is determined by 
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the relationship between s and r . Taking into account the requirement for a constant 
vesicle volume, the following relationship for s is obtained: 

s 2 

s = r - -2- . (3) 
5r 

The correction for a constant volume ([3]) contains a higher term in the powers of S2 and 
thus does not affect in the lowest term either the bending energy or the energy due to 
the electric field, yielding s = r an adequate approximation. 

For small deformations, quadrupolar deformation only induces small perturbative 
changes in the vesicle bending energy and the energy due to the electric field. 

1 9 2 Gbend 



G ! bend('S2) ~ Gbcnd(S2 — 0) + 



2 ds 2 , 

G fle id(s 2 ) « G fie i d (s 2 = 0) + 



2 

<9G n cld 



ds? 



si (4) 

s 2 =0 

(5) 

2=0 

The total energy of the vesicle formally depends on two parameters, ro and s 2 . The 
constraint requiring a constant vesicle volume however eliminates one degree of freedom, 
yielding (j3HHj). It is also worth noting that the fact that neither prolate (s 2 > 0) nor 
oblate shapes (s 2 < 0) have a bending energy lower than those of a sphere means that 
the expansion for the bending energy contains no linear term. 

Equilibrium vesicle deformation, expressed in terms of s 2 , can then be calculated 
by minimizing the total free energy over s 2 : 



—j — (Cbend + Afield ) ~ 7T^ 

dso os z 



S2 + 



ds 



2 



= 0. (6) 

s 2 =0 



'2 

Equation gives the extent of deformation (s 2 ) at given conditions. 

It is now our task to write an expression for the vesicle bending energy (j3J), which 
can be written as [5]: 

Gbcnd = ~k c j> (ci + c 2 ) 2 dA . (7) 

Here, k c is the bending elastic modulus of the membrane, while c\ and c 2 are the principal 
curvatures of the membrane. The spontaneous curvature Co has been omitted, because 
it vanishes for a bilayer composed of two equal layers. The integration is conducted over 
the total membrane area A of a quadrupolarly deformed vesicle. Up to quadratic order 
terms in s 2 , the total bending energy of a nearly spherical vesicle can be written as [5]: 

GbeBd = S7Tk c + ^k c ^ . (8) 

This can be readily interpreted as the bending energy of a sphere plus an addition due 
to the quadrupolar deformation. 

The free energy term arising from the electric field (j3j) is calculated as the work 
done by the forces of the electric field during the deformation of the vesicle. In order 
to evaluate it, we first have to calculate the electric field around a vesicle (figure [2]). 
Gauss' law V ■ D = together with the requirement for an irrotational electric field 
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Figure 2. A vesicle is modelled in spherical coordinates (r, 9) as a spherical shell 
with the external radius r out = ro + d/2 and the internal radius rj n = rg — d/2 
(d <C rin being the membrane thickness) exposed to an external electric field E. 
The conductivity and permittivity of the aqueous solution inside and outside the 
membrane are denoted by a w and e w , and the membrane conductivity is denoted by 
a m . Membrane permittivity is in general anisotropic, hence e\\, ej_; their orientation is 
defined in the text (equation 1 The three media, vesicle exterior, vesicle membrane 
and vesicle interior are denoted by (1), (2), and (3), respectively. 



V x E = 0, which stems from Faraday's law for electromagnetic induction, yields the 
Laplace equation for the electric potential U in a medium with homogenous dielectric 
permittivity [23] - in our case the aqueous medium in the vesicle interior and the external 
aqueous medium (respectively denoted by (3) and (1) in figure |2]). In the anisotropic 
case, however, Gauss' law does not lead to the Laplace equation for the electric potential 
U, V • D = V • eE = —V • (eVLQ. Here, D is the displacement field, E is the electric 
field and e is the permittivity tensor. 

The boundary conditions on both the outer and the inner membrane-water 
boundary impose VxE = 0, thus requiring the continuity of the tangential component of 
the electric field, and that the total surface charge density, including both free charge and 
displacement charge, must vanish. In the spherical geometry, this yields the equations 



for the (1-2) interface (r = r out , figure [2]) and similarly for the (2-3) interface (r = r m ). 
Apart from (|9 HT0|) . the system is constrained by two additional conditions: one requiring 
that the electric field far from the vesicle is unperturbed, and the other requiring that 
the electric field is finite inside the vesicle. 

In (19][T0~I) we introduced the notation for the dielectric permittivities and electric 
conductivities in the three media: e* l > = = e w , = = a w , amd = o m . The 




(9) 
(10) 
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membrane permittivity is treated as anisotropic. Due to its structure, electric properties 
in the direction along the long axis of phospholipid molecules (or normal to the vesicle 
membrane) are expected to be markedly different from the properties in the direction 
perpendicular to it. The phospholipid molecule is not axially symmetric per se, but due 
to the unordered liquid-like structure of the bilayer we can approximate it as such on the 
timescale of interest to our problem. The two directions perpendicular to the normal 
on the membrane can thus be treated as being equal. The component along the normal 
to the membrane is denoted by eii and the permittivity in the direction perpendicular 
to it by e±. Local permittivity in a given point on the membrane is thus equal to en in 
the normal direction and e± in the tangential direction. Permittivity tensor is locally 
defined as: 










e_L 







e± 



For a spherical shell, local coordinates conveniently coincide with the spherical 
coordinates. Gauss' law V ■ (eV£7) = for the treated geometry is most conveniently 
written in spherical coordinates, 



1 d 



II r 2 Q r \ Q r 







'■ sin 9 39 



sm9 le 







(12) 



In spherical coordinates, the equation can be separated into the radial and the angular 
parts, U(r,9) = R(r)Q(9) (it has been already taken into account that the problem is 
axially symmetrical and thus independent of </>), yielding two separate equations for the 
radial and the angular part. For the angular part, a Legendre differential equation is 
obtained. The equation obtained for the radial part differs from the spherical Bessel 
equation by involving a factor e_i_/eii in the second term: 



1 d 

r 2 dr 



2 dR\ 
dr J 



e±l(l + l] 



R = 



(13) 



It is worth noting that ( fT3j) represents a special case of Heun's equation, which arises 
from solving the case of a general ellipsoid [13J. Equation (fl3l) is solved by a linear 
combination 



Rl ( r ) = c ir |(-i-V 1+4 ^/ e ii+ 4p ^/ e ii) + c^K-i+vA+^-lAii+^Aii) . 

The symmetry requirements of our problem limit us to the case 1 = 1. Thus, the ansatz 
for the electric potential that satisfies Gauss' law for the anisotropic case is 



= i[( a (2) r (a-l)/2 + 6 (2) 7 .(-«-l)/2) cc .sfle- 



-iwt 



+ C.C. 



(14) 



A shorthand notation a = ^Jl + 8e_i_/e|| has been introduced along the way. The 
isotropic case corresponds to a = 3, which simplifies (j!4p to a known form 



a {k) r + 



bV>)\ 



J cos 9 e 



—iu)t 



c.c. 



k= 1,3 



(15) 
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The coefficients a^ k \ b^; k = 1,2,3, can in general be complex to allow for a phase 
shift, and are determined from the boundary conditions. 

The boundary conditions for the unperturbed field far away from the vesicle and 
the finite field inside the vesicle immediately yield two coefficients: 

a« = -E Q (16) 

6 (3) = (17) 

The remaining four coefficients are determined by the four equations specifying 
boundary conditions (pi [10]) : 

(2) = -QE r^- a)/2 (2 + (3 + a(3) 

° " 8 + (2 + 6a)(3 + (a 2 - 1)(3 2 - [8 + (2 - 6a)(3 + (a 2 - 1)(3 2 ] 7 * 1 J 

G (3) = -12E ^ 7 (-3)/ 2 

8 + (2 + 6a) /3 + {a 2 - 1)(3 2 - [8 + (2 - 6a)/? + {a 2 - l)f3 2 } 7 * 1 ' 

,(i) = E r 3 [(a-l)(3-2}(2 + f3 + af3)(l-r) , , 

' " 8 + (2 + 6a)(3 + {a 2 - l)f3 2 - [8 + (2 - 6a) + {a 2 - 1)(3 2 } j a 1 ' 

h m = -6E r^l\a-l)P-2]T , . 

" 8 + (2 + 6a)(3 + (a 2 - 1)(3 2 - [8 + (2 - 6a)f3 + (a 2 - l)[3 2 ] T 
Consistently with |8J, two more shorthand notations have been introduced: 

a m — ioje\\ , . 

f3 = — . \ (19) 
a w — iue w 

7 = — ■ (20) 

r out 

The surface density of the force exerted on the boundary of dielectrics by the electric 
field is equal to the scalar product of the Maxwell stress tensor and the normal vector 
to the membrane, 

fout = (T (1) -T( 2 >)e r , (21 fl ) 

f in = (T (2) -T (3) )e r . (215) 

The force vanishes in a homogeneous medium, but can in general be non-zero on the 
boundaries of media with different electrical properties. The Maxwell stress tensor is 
defined as 

T = D®E-i(D-E)J, (22) 

where J denotes the identity matrix. 

Unlike in the case of the bending energy term in the total free energy (jHJ), for 
which an analytical expression was obtained, an approach where energy difference is 
computed is employed here. £Gfi e id denotes a small change in the energy due to the 
electric field, when a sphere (s 2 = 0) is perturbed by a small deformation change 5s 2 - 
This energy difference is calculated as the work done by the forces of electric field during 
the displacement of membrane elements Sr (equations 12 q||2 b\\ . integrated over the entire 
membrane area [8]: 

5G fic i d = - /(font ■ 5r) dA out - /(f in • St) dA m . (23) 
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The integration is conducted over a sphere, which is consistent with the limit of small 
deformations (s 2 "C r ). 

Substituting the coefficients ( I18al - ll8dl) into ( 123]) yields a lengthy expression for 
^Cficid, which will not be reproduced here. After substituting the expresson ([19]) for (3 
into it, one obtains an expression of the form (A + Buj 2 + Coo 4 ) /(D + Euj 2 + Fuj a ) , which 
can be written as a sum of two dispersion terms: 

rfCfidd SG field 6ir p2 2 | f £l , ^2 



rf S2 ' S 2 5 + 1 + ^2 + 1 + U 2 T 2 ) » ( 24 ) 

where the coefficients £oo, £i, £ 2 , Ti and T2 are rather lengthy expressions involving 
five different material constants: en, e±, e w , o~ m , o~ w , the membrane thickness d and 
the vesicle radius r . It turns out, however, that the actual number of independent 
parameters is lower. By introducing a dimensionless frequency u/(a w /e w ), one can 
reduce the frequency dependence of (|24p to only four parameters. Here, we have chosen 
them to be 7, e\\/e w , Ae/e w and a m /a w . 

Equilibrium vesicle deformation can finally be calculated by minimizing the total 
free energy (jBJ) over s 2 , yielding 

. _ 1 4e w E 2 f 6 6 \ , , 

S2 " 16^7" v e °° iT^vf + I+^fJ ' (25) 

where rgr 2 ut has been replaced by r A , consistent with the expansion in E 2 and retaining 
the terms with the lowest energy. As we can see, the dependence of vesicle deformation 
S2 on the circular frequency uj = 2-rv contains two dispersion terms of the Maxwell- 
Wagner type, which arise due to the interfacial polarization and not due to intrinsic 
dispersion. It is also worth emphasizing that the only approximation used in deriving 
the expression ([25]) is that of a small deformation (s 2 <C r ). 

The coefficients £oo, £1, £ 2 , T i an d t 2 figuring in (|2"5T) can be simplified using 
expansion in (1 — 7), as (1 — 7) < 10~ 3 (the vesicle radius is in the range 1-100 \xm, 
while the thickness of the phosphilipid membrane is approximately 4 nm). Expansion 
up to the first order in (1 — 7) yields: 

2[-2e\\{e w - ey) 2 + (e w + 2 e,|)e w Ae + 2e\\{Ae) 
9e|| 2 e L 

-Ae\\a m a w (e w - e») +2a w {e\\a w + e w a m + 2e\\a m )Ae f a m £11 , 

{l= v< (27) 

8(eii + Ae) f a n 



^00 = ^ (1 - 7) , (26) 



6 =7V^^-^ (1 - 7) " (28) 

n =-SL, (29) 



0~m 

e w 



r 2 = ^ . (30) 

We have introduced dielectric anisotropy Ae = e± — en. The characteristic times ri 
and r 2 correspond to the trans-membrane and the trans-vesicle-interior relaxations, 
respectively. It is worth noting that neither T\ nor r 2 is affected by the anisotropy 
in this first-order approximation. 
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3. Results 

This section consists of four points. First, we plot the dependence of deformation on 
frequency and comment on the effects of different material constants. Next, we look for 
the conditions where an oblate deformation (s2 < 0) exists in the high-frequency limit. 
A similar search is performed for the conditions where a prolate deformation (s2 > 0) 
exists in the low-frequency limit. Finally, we plot and analyse the dependence of the 
prolate-to-oblate transition frequency on material constants. 

Dependence of deformation on frequency. Figure [3] shows the vesicle deformation as a 
function of the circular frequency uj of the applied external electric field. Even though 
the expression for deformation ( 1251) contains two dispersion terms, only one of them 
is prominent on the diagram. This is because at values of the parameters used for 
evaluation (ey <C e w , o m <C cr w ), the magnitude of the second dispersion term is about 
two orders of magnitude smaller, £ 2 "C £i- This agrees with the expressions (I27II28I) 
where expansion in (1 — 7) is used: £ 2 oc cr m /a w , while £1 oc (e w /e\\)(a w /a m ). 

Figure [lb provides a comparision of the results of the model with experimental data. 
The parameters used to calculate the dotted line in figure [lb were k c = 0.9 x 10~ 19 J, 
E = 200 V/cm, a w = 10" 4 ^m" 1 , a m = 0.38 x 10~ 5 fl^nr 1 , e w = 80e , e\\ = 2e , 
Ae = 5.2e . While the results of the model reproduce a general trend, it is clear that 
its agreement with the experimental data is only qualitative. This can be attributed to 
the perturbative nature of the model, which is only valid for small deviations from a 
sphere only. This limitation is somewhat conflicting with the experimental requirements, 
where significant deviations from a sphere are clearly preferred for giving a more accurate 
readout. 

The frequency-dependent prolate-to-oblate transition does not reflect a change of 
sign in the forces of the electric field - both in the high- and the low-frequency limit, the 
forces induced by the electric field pull the vesicle apart at poles and compress it around 
the equator, thus favouring a prolate deformation. A lower energy for an oblate shape 
at high frequencies stems from the shear term, fg 5rg. In the high-frequency limit, fg is 
zero for the isotropic case, and increases in magnitude with the increasing anisotropy 
Ae. 

High-frequency limit. High-frequency behaviour depends only on the dielectric 
permittivities of the membrane and the aqueous solution, and is completely determined 
by coefficient ^ as can be seen from f[25l) : if £oo > 0, the deformation is prolate, and 
vice versa. The sign of the exact expression for the coefficient £oo (without employing an 
expansion in (1 — 7)) depends on two parameters; here, they were chosen to be e\\/e w and 
Ae/e w . It is worth noting that while the expression ^ depends on 7 as well, varying 
7 can not change its sign. Furthermore, for all values of vesicle radius which are of 
experimental interest, 7^1 holds, so 7 does not significantly affect vesicle behaviour 
in our model, therefore we did not put much emphasis on it. Values of e\\/e w and Ae/e w 
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Figure 3. The theoretical dependence of vesicle deformation, expressed as either 
s 2/r"o (left axis) or c/a (right axis), on the circular frequency u> of the applied external 
electric field. The curves a, b and c correspond to the membrane conductivities of 
0.5 x 10~ 5 , 1.0 x 10~ 5 , and 2.0 x 10~ 5 fi -1 m _1 , respectively. Other parameters used 
in the calculation were r = 10 urn, 7 = 0.9996, k c = 1.2 x 10~ 19 J, E = 100 V/cm, 
<j w = 10~ 4 n^m -1 , e w = 80e , e» = 2e , Ae = 6e . 




for which ^ = holds can only be computed numerically. In figure HI the transition 
boundary between prolate and oblate shape in the high frequency limit on the (e\\/e w , 
Ae/e w ) plane is plotted. 

An insight into the high-frequency behaviour can also be obtained from the 
expanded expression ( 1261) . If in addition to the expansion in (1 — 7) an expansion 
in eu/e w is employed (it is estimated en/e w ~ 1/40), one obtains the following expression 
for the requirement for an oblate deformation in the high-frequency limit: 

Ae > 2e,| . (31) 

It can also be seen that in the isotropic case (Ae = 0), coefficient ^ attains the 
form (an expansion up to the first order in (1 — 7) is used): 

eoo = 4(e r" £|l)2 (1-7), (32) 
yeye^ 

which is always positive. Thus the deformation of a vesicle with isotropic membrane 
permittivity is always prolate in the high-frequency limit, which is consistent with the 
findings of Winterhalter and Helfrich [8] . 

Low-frequency limit. The condition for prolate deformation in the low-frequency limit 
is obtained from lim^ s 2 , yielding the criterion £00 + £1 + £2 = for the threshold 
between prolate and oblate shapes. Unlike the condition for the high-frequency limit, it 
depends on the dielectric permittivities as well as the conductivities of membrane and 
the aqueous medium. The sign of the exact expression ^00 + ^1 + £2 depends on three 
parameters; here, they have been chosen to be e\\/e w , Ae/e w , and a m /a w . The comments 
on the dependence of the expression on 7 in the high-frequency limit apply here as well. 
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Figure 4. The boundary between the prolate and the oblate shape in the high- 
frequency limit as a function of the ratio between the permittivity of the membrane 
(ey) and the permittivity of the aqueous solution (e w ). Deformation is oblate in the 
area above the curve and prolate in the area below the curve. For eu <C e w , dependence 
can be approximated with a straight line Ae/e w = 2e\\/e w . The curve corresponds to 
a vesicle radius of 10 \xm, yielding 7 = 0.9996. 




Figure 5. The boundary between the prolate and the oblate shape in the low- 
frequency limit as a function of the ratio between permittivity of the membrane ey 
and permittivity of water e w on the one hand and the ratio of membrane conductivity 
<7 TO and the conductivity of the aqueous medium a w on the other hand. Deformation is 
prolate in the area above the curve valid for a particular value of dielectric anisotropy 
and oblate in the area below it. The delineation lines correspond to different values of 
dielectric anisotropy (from bottom to top): Ae = 0, 0.0144, 0.01592, 0.016, 0.018, 0.1, 
3, and 30 eo- The curves are calculated for 7 — 0.9996. 

Its zero can only be computed numerically. In figure O the threshold between prolate 
and oblate shape in the low frequency on the (e\\/e w , a m /a w ) plane is plotted for several 
different values of Ae. The curves are plotted for 7 = 0.9996, which corresponds to a 
vesicle radius of 10 |xm. 

It is worth noting that for small values of anisotropy (Ae « 0), curves intersect 
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the abscissa. For Ae = 0, the intersection point is approximately 1/2 — (1 — 
7)/3, meaning that for ratios e»/e w lower than this value (this corresponds to all 
experimentally achievable cases), the deformation of a vesicle with a membrane of 
negligible conductivity in the low-frequency limit is prolate, while for the cases e\\/e w > 
1/2 — (1 — 7) /3, the deformation in the low- frequency limit for the same vesicle is oblate. 
Increasing dielectric anisotropy lowers the threshold value, and at Ae/e w « 0.0002, the 
deformation of a vesicle with a non-conductive membrane is always oblate in the low- 
frequency limit. This is consistent with the findings of Winterhalter and Helfrich [8], 
who have predicted prolate deformation in the low-frequency limit for a vesicle with no 
dielectric anisotropy and very small membrane conductivity o m jo w = 10~ 10 . 

Again, additional insight can be obtained when one employs expansions in (1 — 7) 
and e\\/e w , yielding the following requirement for prolate deformation in the low- 
frequency limit: 

->t(^-^)' (33 > 

For e_i_ = 3e||, which also fulfills the requirement ( 13TT) obtained for oblate deformation in 
the high-frequency limit, fl33l) simplifies to 

eii 

Cm ~ — o- w . (34) 

In f[3~4"|) . it has been taken into account that eii,ej_ <C e^,. Prolate deformation in the 
low-frequency limit for a vesicle with nonzero membrane anisotropy is only possible 
when the ratio of membrane conductivity and the conductivity of the aqueous solution 
exceeds the ratio of permittivities of membrane and the aqueous solution. 

Another special case is the isotropic case (en = e± -C e^). Here, the expression on 
the right-hand side (I3"3"j) is negative, and so the inequality is always fulfilled, meaning 
that the deformation of a vesicle with no anisotropy in dielectric permittivity is always 
prolate in the low- frequency limit. 



Transition frequency. It is also of interest to plot the transition frequency between the 
prolate and the oblate shape as a function of membrane conductivity. This corresponds 
to finding the zeroes of the following equation for different values of parameters: 



6 6 



«» + r^? + T^ = - (35) 

Parameter space can be reduced by using a dimensionless frequency, as defined above. 
One can calculate the threshold dimensionless frequency numerically. Its dependence on 
the ratio between membrane conductivity and the conductivy of the aqueous medium 
(a m /a w ) for different values of dielectric anisotropy is plotted in figure El It can be seen 
that for values of a m /a w smaller than a certain threshold value - approximately given 
by ( 1341) - one does not obtain a prolate-to-oblate transition at all. 
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Figure 6. The dependence of the transition frequency between the prolate and 
the oblate shape on the value of the ratio between membrane conductivity and the 
conductivy of the aqueous medium (u m /a w ), plotted for three different values of 
dielectric anisotropy (Ae = 4,5,6 eo). The dependence is calculated for e\\ = 2e , 
7 = 0.9996. On the left axis, the scale is plotted in terms of dimensionless frequency 
Wt/(cr w /e w ); on the right axis, the scale is plotted in terms of circular frequency uit for 
a w = 1CP 4 r2 _1 m _1 and e w = 80eo- 



4. Discussion 

The treatment presented in this paper largely follows the outline by Winterhalter and 
Helfrich [5]. It is important to note, however, that substituting isotropic values for the 
coefficients £oo, £i, £2, r i an d T2 into expression (1231) does not lead to the same result as 
equation (21) in their article. The expansions in factors (3 and 7 used by Winterhalter 
and Helfrich show non-analytical behaviour, where consecutive expansions in 7 and 
j3 yield a result different from the one obtained by consecutive expansions conducted 
in the opposite order, i.e., an expansion in (3, followed by an expansion in 7. Using 
Winterhalter and Helfrich's values of material constants, one can see that (3 > 7 holds 
in the high-frequency limit, while in the low-frequency limit (3 < 7 holds. Due to this, 
the frequency dependence of the simplified expression they obtained by series expansion 
differs qualitatively from the exact expression from which the simplified expression was 
derived - while the deformation described by the exact expression decreases with the 
increasing frequency, the simplified expression increases. In this paper, we avoided using 
series expasion in j3 altogether, expressing it explicitely with the frequency uj and the 
material constants independent of frequency. The expression ( 1251) may provide a better 
approximation to the exact expression than the one used by Winterhalter and Helfrich, 
as it retains the frequency dependence of the exact expression, while the expansion they 
used employs an expansion in f3 and is thus actually limited to the low-frequency range 
only. 

In the previous section, we have developed criteria within the framework of our 
model, under which the vesicle shape is prolate in the low-frequency limit and oblate in 
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the high-frequency limit. The criteria we obtained however opened a couple of questions 
about their physical relevance: the tangential component of permittivity is predicted to 
exceed its normal component threefold, the membrane conductivity is predicted to be 
higher than the earlier estimates [8]. 

(i) Membrane permittivity. Within the framework of the presented model, a criterion 
for an oblate shape in the high-frequency limit was obtained, requiring that the 
tangential component of the dielectric permittivity tensor of the phospholipid 
membrane exceeds its radial component by approximately threefold: e± ^ 3eii. This 
is in qualitative agreement with an abundance of experimental evidence claiming 
that the zwitterionic phosphatidylcholine headgroup is oriented almost parallel to 
the surface of a hydrated bilayer and its movement is essentially confined to a 
plane normal to the membrane, employing techniques such as X-ray diffraction, 
neutron diffraction, 2 H nuclear magnetic resonance (NMR) and 31 P NMR as well 
as molecular simulation studies (see [21] and the references therein). 

(ii) Membrane conductivity. Another criterion was obtained for a prolate vesicle shape 
in the low- frequency limit, which required that the ratio of membrane conductivity 
and conductivity exceeds the ratio of membrane permittivity and the permittivity of 
the aqueous medium. In practice, this means a relatively high value for membrane 
conductivity: a m ~ 10~ 5 Sm" 1 . This exceeds not only the value for a thin layer of 
oil, a m ~ 10~ 14 Sm _1 [8j, but also some other estimates for membrane conductivity 
(a m ~ 10~' Sm -1 ) [llj. The high value of conductivity could be attributed to 
the impurities in the lipid and the transient submicroscopic pores in the membrane 
(step 1 in the accepted five-step description of electropermeabilization, [25J). A 
simple model treating the intact membrane and the pores as resistors connected in 
parallel shows that the total area of pores must account for approximately 1/1000 
of the total membrane area, which corresponds to 1000 pores with an area of 
100 A per 1 (xm 2 of membrane area. This significantly - 10 6 -10 8 x - exceeds 
the density of spontaneous pore formation in the membrane [26] . This latter value 
increases, however, with an the applied electric field [27], even though the electric 
field strength in our setup is approximately 70-300 x lower than its threshold value 
neccessary for the pore expansion into micrometer-sized "macropores" [28J. An 
initial membrane tension present in a vesicle with a relative volume very close to 1 
(the estimate for the vesicle shown in figured] is v = 0.9993 ± 0.0007) also works in 
the same direction. Altogether, a high value of membrane conductivity predicted 
by the presented model remains a poorly understood effect, and awaits further 
investigations before a definitive conclusion can be made. 

Geometry considerations. It has been our intention to build the simplest possible 
model which would explain the observed frequency-dependent prolate-to-oblate shape 
transition for small deformations and yield a clear physical picture of the phenomenon. 
The validity of the model decreases once we depart from small deformation. A possible 
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extension of the model lies in a more appropriate description of the vesicle geometry, i.e. 
spheroidal instead of spherical. An effort in the suggested direction might benefit from 
the treatment of a spheroid vesicle with a permeable membrane [29], as well as from 
the calculations of transmembrane voltage for the case of zero membrane conductance 
[301 131] and the solution for the electric potential for the general anisotropic case [T3] . 

5. Conclusions 

The proposed model provides a theoretical explanation for the observed prolate-to- 
oblate transition of phosholipid vesicle shape, when the frequency of the applied AC 
electric field increases, in the case when the electrical properties of the aqueous solution 
inside the vesicle do not differ from the properties of the aqueous solution outside the 
vesicle. Vesicle deformation at varying frequency was calculated numerically. Using an 
expansion into Taylor series we also derived a simplified expression for the deformation. 
The expression retains the frequency dependence of the exact expression and may 
provide a better substitute for the expression that Winterhalter and Helfrich obtained 
by series expansion, which was found to be valid only in the limit of low frequencies. 
The model with the anisotropic membrane permittivity imposes two constraints on the 
values of material constants: tangential component of dielectric permittivity tensor of 
the phospholipid membrane must exceed its radial component by approximately a factor 
of 3; and membrane conductivity has to be relatively high, approximately one tenth of 
the conductivity of the external aqueous medium. Both constraints seem to be justified 
to a certain degree by the claims found in the literature. 
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